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ABSTRACT 

The Lya absorption spectrum of the highest redshift quasars indicates that they are surrounded by 
giant H II regions, a few Mpc in size. The neutral gas around these H II regions should emit 21cm 
radiation in excess of the Cosmic Microwave Background, and enable future radio telescopes to measure 
the transverse extent of these H II regions. At early times, the H II regions expand with a relativistic 
speed. Consequently, their measured sizes along the line-of-sight (via Lya absorption) and transverse to 
it (via 21 cm emission) should have different observed values due to relativistic time-delay. We show that 
the combined measurement of these sizes would directly constrain the neutral fraction of the surrounding 
intergalactic medium (IGM) as well as the quasar lifetime. Based on current number counts of luminous 
quasars at z > 6, an instrument like LOFAR should detect > 2 redshifted 21cm shells per field (with 
a radius of 11°) around active quasars as bright as those already discovered by SDSS, and > 200 relic 
shells of inactive quasars per field. We show that Lya photons from the quasar are unable to heat the 
IGM or to couple the spin and kinetic temperatures of atomic hydrogen beyond the edge of the H II 
region. The detection of the IGM in 21cm emission around high redshift quasars would therefore gauge 
the presence of a cosmic Lya background during the reionization epoch. 

Subject headings: cosmology: theory - intergalactic medium ~ quasars: general - radio lines: general 



1. INTRODUCTION 

The recent discovery of quasars at z > 6 (Fan et al. 2001; 
Fan et al. 2003) and the subsequent observation of two 
Gunn-Peterson (1965) troughs (Djorgovski et al. 2003; 
White et al. 2003) established a new era in studies of 
' the end of the reionization epoch. The absorption of Lya 
, photons from these quasars probes the ionization state of 
the hydrogen in the intergalactic medium (IGM) near the 
quasar redshift and has revealed giant ionized regions (Cen 
& Haiman 2000; Madau & Rees 2000). Comparison of the 
radii of these regions with the ionizing flux implies quasar 
ages of < 10^ years, and an ionization front that is ex- 
- panding with a relativistic speed (Wyithe & Loeb 2004). 

The generation of planned low frequency instruments 
, (such as LOFAB? or the SKA'^) will open a new window 
' on the state of hydrogen in the pre-reionization universe 
through the emission or absorption of 21cm photons rela- 
tive to the Cosmic Microwave Background (CMB) by neu- 
tral hydrogen in the IGM (Scott & Rees 1990; Madau, 
Meiksin & Rees 1997; Gnedin & Ostriker 1997; Shaver et 
al. 1999; Tozzi et al. 2000; Ihev et al. 2002, 2003; Ciardi & 
Madau 2003; Furlanetto, Sokasian & Hcrnquist 2003; Zal- 
darriaga, Furlanetto, & Hcrnquist 2003; Gnedin & Shaver 
2003; Loeb & Zaldarriaga 2003). The giant H II regions 
surrounding the highest redshift quasars should be visi- 
ble via rings of redshifted 21cm emission from warm but 
neutral gas beyond the ionization front. 

The relativistic expansion of the H II region implies that 
the radii measured along and transverse to the line-of-sight 
using the Lya and 21cm techniques respectively, probe dif- 
ferent epochs due to relativistic time delay. In this paper 
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we discuss what will be learned about the state of the 
IGM surrounding the highest redshift quasars by combin- 
ing Lya and 21cm measurements. In § 2 we discuss the 
measurements of radii of the H II region in Lya absorp- 
tion and redshifted 21cm emission. The relation between 
these radii and the neutral fraction of the IGM are then 
discussed in § 3. Some examples of the brightness temper- 
ature profiles around high redshift quasars are presented 
in § 4. In § 5 we discuss the observation of fossil H II 
regions and their utility as a probe of the radiation back- 
ground at high redshift. Finally in § 6 we discuss the 
utility of quasar Stromgren spheres in the calibration of 
the global signature of reionization. § 7 summarizes our 
main results. Throughout the paper we adopt the set of 
cosmological parameters determined by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP, Spergel et al. 2003), 
namely mass density parameters of r2„j = 0.27 in matter, 
rib = 0.044 in baryons, flA = 0.73 in a cosmological con- 
stant, and a Hubble constant of Hq — 71 kms^^ Mpc"^. 

2. OBSERVATIONS OF STROMGREN SPHERES IN 
REDSHIFTED 21 CM EMISSION AND LYa ABSORPTION 

First, we consider the observed signatures of H II regions 
(Stromgren spheres) around the highest redshift quasars 
during their early relativistic expansion. 

2.1. Evolution of the Stromgren Sphere Around a High 
Redshift Quasar 

Neglecting recombinations which do not become impor- 
tant until the expansion becomes sub-relativistic. White et 
al. (2003) derived the evolution of the physical radius of 
the Stromgren Sphere (i?p) as a function of the source age 
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(iage) which can be easily extended to include the possibil- 
ity of a non-zero radius {Rq) at the initial time fage = 0, 



N t 



(1) 



where N is the production rate of ionizing photons by 
the quasar, is the co-moving hydrogen number den- 
sity, xhi is the neutral fraction and c is the speed of light. 
This cubic equation may be solved algebraically (White et 
al. 2003). Alternatively, Rp may be found by integrating 
the differential equation 



dR^ 
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dt ^\N + 47ri?2cxHinO (1 + z f J ' 

with the initial condition i?p(0) = Rq. The latter approach 

(Wyithc & Locb 2004) allows for the inclusion of recom- 
binations. As long as the number of photons emitted per 
time interval At is larger than the number of hydrogens 
in a shell of thickness AR = cAt at Rp, the H I region 
continues to expand at nearly the speed of light. 

Calculations of Rp in this paper assume that the ion- 
ization front is thin. The spectrum averaged mean- 
free path for the ionizing quasar photons is A ~ 
1.5a;^i(l + -z)~^Mpc, which may be compared to the 



bubble radius Rr, 



4.5Mpc, to yield the fractional 



thickness of the ionization front / = (A/i?p) ~ 2 x 

10-3 (i?p/4.5Mpc)"^ x^l [(1 + z)/7.3\~^. The assumption 
of a thin ionization front may be considered reasonable 
provided that .thi ^ 10^^. A low neutral fraction is un- 
likely given the size of the observed H II regions {Rp ~ 
4.5Mpc) around the z > 6.3 quasars (Wyithe & Loeb 
2004). X-ray photons have a longer mean free path but 
are rare and do not contribute significantly to the number 
of ionizations. 

2.2. The Stromgren Sphere Radius Measured Through 

Lya Absorption 

Equations (1) and (2) refer to the Stromgren sphere ra- 
dius at a time t^ge after the source turns on. However 
measurements of the radius through Lya absorption (i?Ly) 
arc made along the line of sight. Lya photons emitted at 
tagc arc not absorbed imtil a time t satisfying the relation 
t = tags + R{t) I c (note that the quasar may no longer be 
active at €). When considering the value of i?Ly measured 
through Lya absorption, the relevant time to consider is 
therefore the age of the quasar at the time when the pho- 
tons being absorbed at the edge of the H II region were 
emitted tage,e = ^age " R{tage) / C- In the absence of recom- 
binations (White et al. 2003), the relation between i?Ly 
and iagc.c is equivalent to the infinite speed of light solution 
for the evolution of a Stromgren sphere radius (Haiman & 
Cen 2000; Madau & Rees 2000), 
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where we add the subscript || to denote the isotropic equiv- 
alent value of N for a measurement done parallel to the 
line-of-sight. 



2.3. The Stromgren radius measured through HI 
emission and absorption 

Madau, Meiksin & Rees (1997), and Tozzi et al. (2000) 
have considered the H I signature of the IGM surrounding 
a high redshift quasar. Outside the H II region, neutral 
hydrogen may be heated by soft X-rays and by scattering 
of Lya photons. The spin temperature of the neutral hy- 
drogen is coupled to the kinetic temperature of the IGM 
through Lya scattering of a background of UV photons, 
and results in 21cm emission (in excess of the CMB flux) 
around the quasar. 

2.3.1. The Lya Scattering rate 

The spin temperature of a warm IGM quickly becomes 

coupled to the CMB radiation, so that it yields no con- 
trast for H I observations (Madau, Meiksin & Rees 1997). 
However scatterings of Lya photons couple the spin tem- 
perature of the warm IGM at radii R > Rp through the 
Wouthuysen-Field eff'ect (Wouthuysen 1952; Field 1958). 
The spin temperature 
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is weighted between the CMB temperature and the IGM 
temperature through the Lya pumping efBciency 
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In equation (4), = 0.068K is the temperature corre- 
sponding to the transition energy, = 2.85 x lO-^^s"^ 
is the spontaneous decay rate of the hypcrfine transition, 
and Pa is the Lya scattering rate per hydrogen. 

The value of Pa may be calculated for a continuum 
UV source embedded in an expanding IGM around a 
Stromgren sphere as follows (see Loeb & Rybicki 1999 for 
more details on the related physics) . In the absence of scat- 
tering, photons that start at a frequency separation Av on 
the blue side of the Lya resonance frequency Va would red- 
shift into resonance at some distance R^es = {Av/va)c/H 
from the source, where H = H{z) is the Hubble param- 
eter at the quasar redshift. At this distance, the Hub- 
ble velocity of the IGM introduces a Doppler shift that is 
equal to the initial blueshift of these photons relative to 
the Lya line center. When a resonant photon enters the 
neutral IGM near i?res, its mean- free-path becomes very 
short because of the large cross-section to Lya scatter- 
ing. As the photon frequency redshifts very close to the 
Lya resonance, the mean-free-path diminishes, and the 
frequency shift per scattering (which corresponds to the 
Hubble velocity increment across a mcan-free-path) goes 
to zero. Consequently, the photon spends a very long time 
around Rres- At this stage, the finite temperature of the 
gas becomes important. The photons are able to diffuse 
out of resonance primarily due to the Doppler shifts they 
encounter by the thermal speed of the atoms at -Rros- 

Each scattering near R-^cs leads to a fractional frequency 
shift by the thermal Doppler effect of Avti^/va = Vth/c, 
where Vth — (SfcTiGM/wp)^/^ is the thermal speed of the 
atoms and Va is the resonant Lya frequency. The thermal 
shifts are of random sign and so the photons go through 
a random walk in frequency. Consider the photons just 
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around the resonant frequency at i?ics- These photons 
change their frequency by ~ Afth in cacli scattering, and 
have a mean free path of Ath = l/[n^(l + z)^crth]- Here 
dth is the Lya cross-section at the center of the thermally 
broadened Lya line (Rybicki & Lightman 1979) 

5.87x10— cm^f , (5) 



where e and rUe are the electron charge and mass, and 
/21 is the oscillator strength for the Lya transition. The 
mean intensity of nearly resonant photons at Rj-cs is in- 
creased by a factor of ~ ^Ath relative to optically-thin 
conditions, where I is the so-called Sobolev length, defined 
as the distance over which the bulk velocity gradient of 
the medium equals the thermal speed of an atom, 

~ dv/dR ~ H' 
and V = HR is the Hubble velocity. The optical-depth 
factor Tth = ll\h. reflects the fact that the diffusion speed 
of photons across the resonance region of width ^ Hs a 
factor Ath/^ times smaller than the speed of light. 

In a steady state, the mean intensity at all radii should 
therefore be enhanced relative to the optically-thin case 
by a factor of ~ Tth- However, a steady state will not be 
established immediately when light from the quasar first 
penetrates the neutral IGM because of the slow diffusion 
speed of the resonant photons. At early times, scattering 
will produce an exponential reduction in resonant photon 
flux by a factor oc e"^"". Even though photons may be 
sufficiently blue that they do not redshift into the reso- 
nance until Rrcs, they would still be absorbed by the blue 
wing of the Lya resonance (where the cross-section is still 
considerable) at radii smaller than i?res. 

Let us therefore consider photons of arbitrary frequency 
on the blue side of the Lya resonance and calculate their 
mean intensity at -Rrcs. We define the mean intensity as 
the angular average of the specific (directional) intensity 

I{v, ij,), namely J{u) — i J^^ dfj,I{i', fi) where /x = cos^ is 
the direction relative to the radius vector at i?res (Loeb & 
Rybicki 1999). The Lya optical depth for photons with 
frequency v a,t a, radius i?res is 

t(i/) = [ dR'aa 
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where (Ta{v) = (Tth exp{ — (i/- 



c 



R')]]nUl+zr 



(7) 

^ / Au^^} is the thermally- 



broadened cross-section (Rybicki & Lightman 1979), and 
Ai?pi = i?rcs ^ Rp{t ^ At) is the path- length traversed 
by a photon through the neutral IGM after it had crossed 
the boundary of the H II region at a time t — At. In the 
strong scattering regime {t{i') » 1), the mean intensity is 
a factor of r(i/) larger than the specific intensity along the 
radial direction, because the radial diffusion speed of the 
photons is ^ c/r(z/). The initial mean intensity of nearly 
resonant photons at radius i?rcs is therefore 
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where 
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is the redshifted source spectrum in erg s~^ Hz~^. 
The scattering rate Pa may then be calculated from the 
mean intensity as 



dv- (Ja[y). 



(10) 



The directional intensity of resonant Lya photons starts 
high near the boundary of the H II region at Rp , but is ini- 
tially suppressed at larger radii. Resonant photons will dif- 
fuse to radii larger than Rp at a speed ~ c/rth- In the case 
of a stationary H II front, the high intensity region expands 
with time, extending a distance (c/rth)At beyond the edge 
of the H II region after a time At, and absorption will 
truncate the intensity at radii larger than Rp + {c/rth) At. 
For the quasars under consideration the expansion speed 
of the H II region (> 0.1c) is well in excess of the diffu- 
sion speed^ of resonant Lya photons (c/rth ~ lO^^c). On 
the timescales of interest here, there is a deficit of photons 
near the line center as described by equation (8). 

The critical thermalization rate above which Lya scat- 
tering couples Ts to TiGM is given by (Madau, Meiksin & 
Rees 1997) 
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We find that Pa only exceeds the thermalization rate 

Pth within a negligibly thin shell around the H II region 
(< 10~^Mpc). As a result, detection of redshifted 21cm 
emission from the H II regions of quasars will require the 
presence of a UV background. Background photons be- 
tween the Lya and the Lyman-limit frequencies can prop- 
agate throtigh the universe and redshift with time to the 
point where they resonate with the Lya line and couple 
Tg to Tk- Such a background is expected to be present at 
z 6.5 since the large electron scattering optical depth 
measured by WMAP (Kogut et al. 2003) implies that the 
IGM had been significantly reionized by that time (requir- 
ing the production of more than one ionizing photon per 
baryon above the Lyman-limit). 

Our result for the scattering rate should be compared 
with the result of Madau, Meiksin & Rees (1997) who es- 
timated 
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In the steady state this expression incorrectly neglects the 
rise in intensity that follows from the small diffusion speed 
of photons near the Lya resonance and hence underesti- 
mates the value of Pa by a factor of ~ nj\(l -I- z)^athVth/H 
which amounts to several orders of magnitude. However, 
in the initial state the expression neglects the absorption of 
resonant photons before they reach a radius R > i?p, and 
so overestimates the value of Pa outside the H II region by 
several orders of magnitude. 

The scattering of Lya photons may also provide a source 
of heating for the gas. As this peak shifts to longer wave- 
lengths, it loses energy to the hydrogen atoms off which 



^ Note that the diflfusion speed of a photon is a function of time as its frequency drifts in and out of resonance with the Lya transition of the 
expanding IGM. 
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the photons scatter (Chen & Miralda-Escude 2003). The 
additional energy density due to the scattering is 



U 



(13) 



The rate at which energy is lost from the hne in an ex- 
panding IGM equals the heating rate of the gas. Hence we 
have 

r„ = ^ = c/i^ = iip„^. (14) 

at V at cCTth 
This heating rate, which is proportional to the Lya scat- 
tering rate, does not contribute to the heating of gas out- 
side the H II region owing to the very low Lya intensity 
there. Instead, the heating is dominated by X-rays which 
we discuss next. 

2.3.2. X-ray Heating Beyond the H II Region 

Gas in the IGM outside the ionized region is subject to 
secondary electron heating resulting from ionizations by 
soft X-rays that propagate into the neutral IGM beyond 
the H II region. The heating rate Fx in erg s~^ Mpc^'' at 
radius R > Rp due to a source with a spectrum e(i/) (in 
erg/sec/Hz) is 

T^{R,t) = £ dvfJ^^LaMxnin%{l + zf 

X exp [-Ai?pic7pi(i/)xHin?^(l + zf] , (15) 
whore (Tpi(i/) is the cross-section for photoionization and 
Vion = 3.29 X lO^^Hz is the Lyman-limit frequency cor- 
responding to the ionization threshold of hydrogen. In 
equation (15) Ai?pi = R — B,p(t — At) is the path- length 
of a photon through the neutral IGM after it crossed the 
boundary of the H II region at time t — At. We have ne- 
glected adiabatic cooling due to cosmological expansion in 
this expression since the quasar lifetime is much shorter 
than the Hubble time. The fraction of photon energy con- 
verted into heat, /x, may be evaluated using the fitting 
formula (ShuU & van Steenberg 1985) 



/x = 0.9971 (l - [1 - (1 - xm) 



0.2663j 1-3163\ 



(16) 



This equation is valid for photon energies hpV > lOOeV. 
At lower energies, the value of /x may be roughly approx- 
imated by multiplying equation (16) with an additional 
factor of exp([100eV— hpv] /lOOeV) (with the constraint 
that /x < 1). The residual ionization fraction from cos- 
mological recombination of ^ 5 x 10^^ yields /x ^ 0.17 
(Madau et al. 1997), which sets the minimum heating rate. 
The heating rate leads to a temperature rise according to 
dTiGM 2 Tx 1 ^^^^ 



dt 
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2.3.3. Redshifted 21cm Observation of Quasar H II 

Regions 

The above discussion has shown that X-rays from the 
quasar are capable of heating the gas outside the H II re- 
gion, and that the spin temperature of this warm gas can 
be coupled to its kinetic temperature in the presence of 
a strong UV (Lya) background. In regions of the IGM 
where the spin temperature is decoupled from the CMB 
temperature, a differential antenna temperature of 
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is observed. In equation (18) we have assumed a uniform 
IGM at the mean density with a neutral fraction xm- The 
luminous high redshift quasars reside in rare high-density 
regions. However on the scale of several Mpc, the overden- 
sity due to infall around these regions is < 10% (Barkana 
& Loeb 2003a). 

The brightness temperature of the redshifted 21cm sig- 
nal is small and subject to severe foreground contamina- 
tion problems. Studies by Di Matteo ct al. (2002), and Oh 
& Mack (2003) have; found that on scales of several arc- 
minutes, randomly distributed point sources would have 
to be removed down to a level of ~l/zJy before bright- 
ness temperature fluctuations of order 10 mK could be 
detected. In the presence of the expected level of clus- 
tering of these sources, even the removal of point sources 
below 1/uJy would leave foreground fluctuations on this 
scale that are several orders of magnitude in excess of the 
redshifted 21cm signal. However, the spectrum of the fore- 
ground contaminants is expected to be smooth (Shaver et 
al. 1999). The width of the H II region corresponds to a 
recessional velocity v < 3000 km s~^ or a redshift interval 
of ~ {1 + z)v/c ~ 0.07. The redshifted 21cm ring will 
therefore produce a bump (or a dip) on top of the smooth 
foreground spectrum at the redshift of the quasar. 

At a frequency v^p ^ 120MHz, an instrument like LO- 
FAR is expected to have a band-pass of Ai^p ~ 4MHz. 
This should be compared with the frequency shift associ- 
ated with the redshift difference between light emitted at 
the quasar and the edge of the H II region along the line 
of sight 
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Since many channels are contained within each band-pass, 
the structure of the H II region can therefore be probed 
along the line-of-sight in frequency space (Madau, Meiksin 
& Rees 1997), allowing determination of the full three di- 
mensional shape of the H II region, including the radius of 
the H II region behind the quasar in addition to the radius 
i?Ly of the H II region in front of the quasar. 

There are two distinct possibilities for the redshifted 
21cm signature of a high redshift H II region. First, the 
region might be expanding into a cold IGM, in which 
case emission will only be seen out to a radius where the 
warming X-rays can reach. Alternatively, the IGM may 
have been pre-heated, either by an X-ray background, by 
gravitationally induced shocks in the IGM (Furlanetto & 
Loeb 2003), or by an early reionization (e.g. Wyithe & 
Loeb 2003a; Cen 2003). The presence of Gunn-Peterson 
troughs in quasar spectra at ^; > 6 implies an IGM that is 
at least partially neutral, while the relatively small sizes of 
the Stromgren spheres around these quasars suggest a neu- 
tral fraction that is of order unity (Wyithe & Loeb 2004). 
On the other hand, observations from WMAP (Kogut et 
al. 2003) suggest an optical depth to electron scattering of 
T ~ 0.17 and therefore significant reionization at z ^ 20. 
If the universe were reionized twice, then the Stromgren 
spheres would be expanding into a partially neutral, but 
warm IGM. 
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Fig. 1. — Schematic illustration of the geometry of Lya and redshifted 21cm observations of a relativistically expanding quasar Stromgren 
sphere (dark regions). The four panels show the evolution at 4 different times. The quasar turns on at to, emitting photons isotropically 
into a pre-existing H II region of radius i?o- Lyo? photons (open arrows) emitted at to reach the edge of the H II region at t\. At that time, 
21cm photons (closed arrows) are emitted by the heated gas perpendicular to the line-of-sight (at the quasar redshift) at i? > Ro and travel 
towards the observer. At tg, the Lya photons emitted at t\ reach the edge of the H II region. The absorption of these photons encodes a 
measurement of Ro,. When these photons are observed, they are accompanied by 21 cm photons (emitted at t\ from the quasars redshift), 
which show a ring of radius Rq. 



3. COMPARISON OF STROMGREN SPHERE RADII 
MEASURED USING LYa AND H I OBSERVATIONS 

The most natural place to search for the redshifted 21cm 

signature of quasar H II regions is around the highest red- 
shift quasars discovered in wide field surveys such as the 
Sloan Digital Sky Survey (SDSS^). Such quasars have H II 
regions that are also probed via their Lya absorption spec- 
tra. In this section we discuss the complimentary nature 
of the two observations with respect to the determination 
of the quasar lifetime, tage, and the neutral fraction of the 
IGM, xm- 

Stromgren spheres expand faster into a partially ion- 
ized IGM. If the IGM has a neutral fraction xhi, then the 
observation of i?Ly through Lya absorption measures the 
combination (fage,e/a^Hi)^^^ (where here we have neglected 
recombinations and assumed Rq = for simplicity). De- 
tection of the redshifted 21cm ring reflects the state of the 
H II region at tage,e- The differential antenna temperature 
(5Tb is proportional to xm due to the small optical depth 
at 21 cm, 

6n cx (1 - e-"™) « THi (X xm. (20) 

Measurement of the contrast in the emission and absorp- 
tion rings, or of the brightness temperature of the emis- 
sion ring, will facilitate a measurement of xhi independent 
from fage,e- Redshifted 21cm observations may therefore 
provide an opportunity to break the degeneracy between 
tage,e and Xm- Unfortunately, the heating rate and tem- 
perature of the IGM are uncertain, and unless Tk ^ Tcmb, 
the value of (5Tb will be sensitive to the unknown value of 
Tk. Below we discuss an alternative geometric method for 
finding .thi that is not sensitive to the precise value of (5Tb. 

The relativistic expansion and finite light travel time 
imply that the 21cm and Lya measurements sample two 
different epochs in the evolution of the H II region. The 
situation is illustrated schematically in figure 1. Suppose 
that the quasar turns on at to, emitting photons isotropi- 
cally into an existing H II region of radius Rq . Lya photons 

^ See http://www.sdss.org/ 



emitted at will reach the edge of the H II region at ti . At 
that time, 21cm photons with (5Tb > are emitted by the 

heated gas at R > Rq and travel towards the observer. For 
simplicity consider 21cm photons that are emitted from 
the heated IGM just beyond the ionizing front and per- 
pendicular to the line-of-sight (i.e. at the quasar redshift). 
The sphere continues to expand at nearly the speed of light 
between ti and t;^. At t^, the Lya photons emitted at ti 
reach the edge of the H II region. The absorption of these 
photons encodes a measurement of R3. The continuum 
and 21cm photons propagate to the observer, who mea- 
sures a radius R3 along the line-of-sight in Lya absorption, 
and a different radius Rq in redshifted 21cm photons. 

In this paper we define i?Ly to be the radius measured 
through Lya absorption along the line-of-sight. This is 
compared with Rm, defined to be the radius of the ring 
measured perpendicular to the line of sight through red- 
shifted 21cm emission. The radii i?Ly and Rm are not 
equal for a spherical region due to finite light travel time 
effects mentioned in the previous paragraph. We note that 
this effect is separate from the supcr-himinal transverse 
expansion of the ring of 21cm emission from the relativis- 
tically expanding portions of the H II region near the line- 
of sight. These super- luminally expanding rings would be 
seen ^ 3000km s^^ blueward of the quasar redshift (for 
Rp - 4Mpc). 

The exact redshift of the quasar may be uncertain by 
up to 1000km s~^ (Richards et al. 2002). However this 
uncertainty corresponds to only 30% of the diameter of the 
H II region (~ 3000km s~^). A point on the edge of the 
H II region that is blueshifted by ^ 1000km s^^ relative to 
the quasar lies in a direction away from the quasar that 
makes an angle of /9 ~ cos~^(1000kms~^/3000kms~^) ~ 
70° with the line-of-sight. As a result, the misidenti- 
fication of the quasar redshift will lead to an error of 
~ i?p(l — sin/3) ^ 0.05i?p in the transverse diameter of 
the ring at the quasar redshift. Moreover, since P > 70°, 
the ring of 21cm emission at the observed quasar redshift 
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Fig. 2. — The neutral fraction as a function of the ratio Kni/^Ly calculated from equation (24). The calculations are for the case where 
iJLy = 4.5Mpc, z = 6.5 and iV|| = lO^'^s-i. We show curves for N^/N^^ = 1.0,0.3 and 0.1, and for Ro = (left hand panel) and Ro = 2 
(right hand panel). 



will not be subjected to superluminal expansion, even if 
the quasar redshift is misidentified by ~ lOOOkms"^. 

More quantitatively, the relativistic expansion and finite 
light travel time lead to several interesting possibilities: 

• The coupling of the spin temperature of neutral hy- 
drogen to the temperature of the warm IGM sur- 
rounding the H II region results in redshifted 21cm 
radiation in excess of the CMB. _Rhi is measured 
transverse to the line-of-sight and so corresponds 
to the radius at the time the observed UV photons 
were emitted. We can directly infer the ionizing 
photon rate (iVy) which drives the evolution of the 
H II region along the line of sight. However, since 
the quasar emission may not be isotropic, we al- 
low the evolution of the H II region transverse to 
the line of sight to be driven by a different ioniz- 
ing photon rate N±. Thus, the ratio iVj_/iV|| is a 
measure of the anisotropy in quasar emission. Since 
high-redshift quasars are brightness selected, we as- 
sume that N±/N^\ < 1. Using equation (3), we find 
the age of the quasar when light crossing the edge of 
the H II region was emitted (iage.e) and substitute 
it into equation (1) for the radius of the H II region 
at that time. This yields the following relationship 
between i?Ly and Rm 



p3 p3 



+ - 



1 



-Rhi 3 



i?3 



HI 



^1 



(21) 



or 



fRmX 

\RLyJ 




Note that equations (21) and (22) arc indepen- 
dent of Rq when the quasar emission is isotropic, 
and that the dependence of Rni/Rhy on iV_L/iV|| 



is rather mild. If the quasar emission is isotropic, 
the transverse radius Rui is always smaller than 
the line-of-sight radius Rhy, particularly during the 
relativistic expansion phase. When the expansion 
becomes sub- relativistic, i?Hi approaches i?Ly 

• The radii measured through Lya absorption and 
21cm emission correspond to different epochs sep- 
arated by Rhy/c. For a spherical H II region we 
obtain the average expansion velocity of the ioniz- 
ing front between the times corresponding to Rm 
and i?Ly 

,dRp ( Rm\ 



(23) 



Note that this velocity is the physical velocity of the 
ionizing front. The apparent transverse velocity of 
portions of the ring might be superluminal. 

• Given a fixed luminosity, the expansion speed at 

a radius J?Ly is proportional to x^^^ , and so the 
neutral fraction can be inferred from the two radii. 
Equation (22) yields 



a;Hi 



R 



'Ly 



c n^(H-z) 



R 



HI 



N± 



Rq 



3. 

N± 



(Rm/RLy] 



\Ri.yJ \Rl.yJ {WI 

V 



R 



Ly 



,3.25Mpcy 




Some examples of the dependence of xm on the 

ratio (Rm/Rhy) given i?Ly = 4.5Mpc and 7V|| = 
lO^^s"^, are shown in figure 2. In the case of 
isotropic quasar emission (A^_L/iV|| = 1), a small 
ratio of (i?Hi/^^Ly) implies relativistic expansion, 
and hence a small neutral fraction unless it!Ly ^ 
3.25Mpc. Values of A''j^/A''|| which are less than 
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Fig. 3. — Plots of the heating rate of the IGM dTiQ^/d-t (in K yr^-'^) due to the X-ray emission by the quasar, as a function of the Stromgren 
radius. The three curves correspond to radii at which the photons have traversed H I path lengths of ARpi = 0.01i?p, 0.1i?p and 0.9i?p, 
where Rp is the physical radius of the H II region. The assumed spectrum is given in equation (25), and the neutral fraction xm was set to 
unity. 



unity lead to larger inferred values of Xfji- Simi- 
larly, a non-zero value for Rq leads to larger values 
of Xm- The assumptions of isotropy and Rq = 
lead to determination of a lower limit for xm- The 
observation of a redshifted 21cm ring with a radius 
greater than 2.5Mpc (0.5Mpc) around the known 
z > 6.3 quasars (with i?Ly ~ 4.5Mpc) would pro- 
vide a limit on the neutral fraction of Xm > 0.3 
{xm > 0.05). 

4. SOME EXAMPLES 

In this section we present examples of the redshifted 
21cm signatures of 2 > 6 quasars that will be observed by 
instruments like LOFAR. The UV continuum blucward of 
Lya is assumed to be that of the median spectrum in the 
sample of HST quasars studied by Telfer, Zheng, Kriss & 
Davidson (2002). They found a power-law €c{i') cx //""euv 
with aEUV = 1-57. This slope may be extended into the 
X-ray regime based on the results of Yuan, Brinkmann, 
Siebert & Voges (1998). We include the Lya emission 
line centered at = 912A, with an equivalent width of 
Wx = 90A(Vandcn Berk ct al. 2001; Tclfcr, Zheng, Kriss 
& Davidsen 2002), and a line width of (j\ = 19A(Vanden 
Berk et al. 2001), providing a UV-X-ray spectrum 



e{v) = 1.3 X lO-^^ergs^^Hz 



1 



■ exp 



[cjv - A„)^ 



(25) 

The normalization of this spectrum is appropriate for a 
2 X 1O^M0 supermassive black hole accreting at its Ed- 
dington rate (Elvis et al. 1994), and it yields an ionizing 
photon rate of TVy = 1.3 x 10^''s~"'^, representative of the 
z > 6.3 SDSS quasars (White et al. 2003). In this section, 
we show examples of the 21cm emission signal that may 
be observed around the known very high redshift quasars, 
with the spectrum described by equation (25). 

The level of redshifted 21cm emission is enhanced 
through heating of the IGM by the quasar X-ray emis- 
sion. The corresponding heating rates dTiGu/dt of the 
above quasar spectrum are plotted as a function of ra- 
dius in figure 3. The three curves correspond to radii 
where photons traversed path lengths through neutral gas 



of Ai?pi = O.OliJp, 0.1i?p and 0.9i?p. The neutral fraction 
was set to unity in these examples. Note that lower neu- 
tral fractions lead to higher values of /x and hence larger 
heating rates. 

If the Stromgren spheres expand into a pre-heated IGM 
containing a UV background that coupled the spin tem- 
perature of the IGM to its kinetic temperature, the H II 
regions should appear as holes in the sky of redshifted 
21cm emission. As an example we first consider an IGM 
that is pre-heated to ~ lOOK at 2; = 6.3 (Tigm ~ STcmb)- 
We follow the evolution of the Stromgren sphere, and at 
each time compute the heating rate of the IGM beyond 
the edge of the H II region using equation (15). The evo- 
lution of the kinetic temperature of the IGM is followed on 
a grid using equation (17). The spin temperature is then 
set to TiGM as appropriate in the presence of a strong Lya 
background. Figure 4 shows four examples of the 21 cm 
emission signature of a Stromgren sphere expanding into 
a warm IGM as a function of radius (top axes) and ap- 
parent angle (bottom axes) for different combinations of 
N±/N\\ = 1.0,0.1 and xm = 1.0,0.1. An initial value of 
Rq = was assumed. The y-axcs arc labeled by the dif- 
ferential antenna temperature ST^ (right axes), and the 
flux SI (left axes) assuming the design goal sensitivity of 
LOFAR 

-2 



61 = (1.2/uJyarcmin ) 



1 



lOmK/ V 7.5 



(26) 



and a top-hat beam with a diameter of 2'. In each case, 
four profiles are shown for observing times correspond- 
ing to quasar ages of tage.e = 0.5,1,2,4 x lO^yr at the 
time when the optical photons observed along the line-of- 
sight were emitted. The lower row includes cases with 
Xm = 0.1. The smaller neutral fraction results in larger 
radii, faster relativistic expansion, and therefore thinner 
emission shells. The 21cm emission signal is weaker in this 
case because ST^ oc .thi- There is excess emission above 
the level of the general IGM in a ring covering the region 
where gas has been heated by the quasar. However the 
brightness temperature is only slightly higher than the rest 
of the IGM. This is due to the temperature dependence of 
ATb which asymptotes to a constant value of ^ 23a;HimK 
for high temperatures. The panels on the right show exam- 
ples with iV_L/iV[| = 0.1. The lower transverse flux results 
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Fig. 4. — Four examples of the evolution of a Stromgrcn sphere expanding into a warm (Tigm = STcmb) IGM as seen in redshifted 21cm 
emission and Lya absorption, assuming different combinations of N±/N^^ = 1.0,0.1 and x-jji = 1.0,0.1. In all cases wc assume Ro = 0. The 
21cm emission signal is plotted as a function of radius (upper horizontal axes) and apparent angle in arc-minutes (lower axes). The y-axcs 
are labeled by the differential antenna temperature 5Tb (right axes) and the flux SI (left axes) given a circular top-hat beam with a diameter 
of 2'. In each case, four profiles are shown for observing times corresponding to quasar ages of tage,e = 0.5, 1, 2,4 x lO^yr when the optical 
photons were emitted. Values of Rj^y at tage.e = 0.5, 1, 2, 4 x 10''^yr are shown with tick marks. 
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Fig. 6. Plots of the radius measured through Lya absorption along the line-of-sight -RLy(tagc) (solid dark lines) as well as the radius of 
the interior boundary of the 21cm emission ring at the quasar redshift (solid light lines). The dashed lines show the speed at which the edge 
of the H II region seen in 21cm (light lines) and Lya absorption (dark lines) is moving when observed at a time corresponding to tage.e- The 
curves correspond to the example shown in the top left panel of figures 4 with Arx/iV|| = 1 and x^i = 1. 



in an expansion that becomes sub-relativistic at earlier 
times. As a result, the emission shell has a smaller radius 
and a larger thickness. 

As our second example wc consider an IGM that is cold 
in neutral regions [T ^ 0.026(1 + z)"^, corresponding to 
adiabatic cooling following the thermal decoupling of the 
baryons from the CMB]. The resulting 21 cm signal is 
shown in figure 5. The heated portions of the IGM beyond 
the edge of the H II region show emission, while the cold 
IGM produces absorption of the CMB at larger radii. The 
brightness temperature of the emission ring obtains values 
of Tb ~ 20a;HimK, smaller in amplitude than the absorp- 
tion signal of Tb ^ — 295xHimK. A cold IGM therefore 
offers the best contrast for observation of the transverse 
extent of the H II region. As in the previous example, the 
emission rings get thicker at late times because the expan- 
sion speed becomes less relativistic. Since the transverse 
expansion becomes sub-relativistic sooner than the line- 
of-sight expansion, the rings become thick at earlier times 
for a lower value of N±/N^^ . 

Figure 6 shows i?Ly(iagc) (solid dark lines) as well as the 
inner radius of the 21cm emission ring (solid light lines) 
for the example presented in the top left panels of figures 4 
and 5 {N±/N^\ = 1, .xhi = 1). The four epochs correspond- 
ing to the differential antenna temperature profiles plotted 
in figures 4 and 5 are indicated by the vertical tick marks. 
Similarly, values of i?Ly at tage.e = 0.5, 1,2,4 x lO^yr are 
indicated by tick marks in the panels of figures 4 and 5. 
The value of i?Ly is always larger than the radius of 21cm 
emission. The dashed lines show the velocity at which the 
edge of the H II region is moving when observed at a time 
corresponding to tage.e- The 21cm ring (light dashed lines) 
is observed at an earlier phase of evolution and so is mov- 
ing at a higher velocity. The edge of the H II region is 
moving at a speed > 0.2c when Rp < 4.5Mpc. 

5. FOSSIL H II REGIONS IN THE HIGH REDSHIFT IGM 

The presence of a Lya background has very interesting 
implications for the observation of fossil H II regions. The 

recombination time within an ionized region of the IGM 
at mean density is ^rccom ^ («b'^h(1 + z)^)^^ ^ 1-5 x 
10^[{1 + z)/?.^]-^ years, where = 2.6 x lO^^-'^ cm'"* s"! 
is the case-B recombination coefhcient at a temperature of 
lO^K. Thus, the recombination time at the mean cosmic 



density exceeds the Hubble time for z <7. Once a quasar 
has turned off, it therefore leaves behind a fossil H II re- 
gion that survives recombination. This H II region will 
therefore remain detectable in redshifted 21cm radiation 
in both scenarios discussed in this paper. Since the bright 
SDSS quasars are likely to be active for only 1% of the 
time (Wyithe & Loeb 2003b), there should be of order 
100 times as many 21cm emission rings around fossil H II 
regions as around active quasars. Furthermore, the hier- 
archical buildup of the supermassive black holes powering 
the quasars would result in an even larger number of small 
fossil H II regions. 

How many H II regions should an instrument like LO- 
FAR detect? It is anticipated that near frequencies of 
Z4,p ^ 200MHz, the virtual core of the planned instrument 
LOFAR will have a beam-size of A6'bcam ~ 23°, a resolu- 
tion of several tens of arc-minutes, and a system tempera- 
ture of Tgys ~ 300K. The instrument will have an effective 
collecting area of ^ ~ lO^m^ at this frequency, leading to 
a resolution of A6'fa ^ 1.22AyV/]4 ~ 22' for the telescope 
in a filled aperture mode. Higher resolutions are obtained 
by spreading the array out over longer baselines, at the 
cost of decreased sensitivity. From the radiometer equa- 
tion (e.g. Burke & Graham-Smith 1997) we can find the 
root-mean-square noise in the brightness temperature as a 
function of resolution on angular scales A6 < A9fa, of 

Ar.,A.) = (0.24„K)(^)^(|^) 

J^Y't^Y"' (27) 

V4MHzy VlOOhr/ ^ ' 

(28) 

Here Afbp is the detector band-pass, and r is the inte- 
gration time. The noise in the brightness temperature is 
constant on scales larger than A9{e,. The root-mean-square 
noise for an observation with a resolution A9 ^ 3' there- 
fore becomes comparable to the expected emission signal 
((5Tb ~ lOmK) for a modest integration time of r ~ lOOhr. 

The co-moving optical luminosity function of bright 
quasars at z ~ 6 in units of Tq^^Gpc^'^ (Fan et al. 2003) is 

$(L) « lO-^'L-'Gpc-' ( loT^) ' ' ^'^^ 
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where P ^ —3 is the logarithmic slope at the bright end. 
Given a band-pass Au^p ~ 4MHz and a central frequency 
Z4,p ^ 200MHz, the number of active quasars per loga- 
rithm of luminosity per field of view is 

where dPV /dz dfl is the co-moving volume per unit red- 
shift per unit solid angle. Equation (30) yields a value of 
N 2 bright ~ 6 quasars per LOFAR field. 

The number of redshifted 21cm rings or holes in the 
emission signature of the IGM per field as determined by 
equation (30) applies to active quasars. However as noted 
above, the 21cm emission may persist after the quasar has 
turned off. This is because in the absence of other heating 
sources, gas heated to temperatures far in excess of the 
CMB will cool adiabatically. The (1 + z)^ dependence of 
the temperature in this regime should ensure that the gas 
remains at a temperature in excess of the CMB by the time 
it is re-heated during the next episode of quasar activity. 
Wo might expect to find significant numbers of fossil H II 
regions, provided that there is continued coupling of Tg 
and TiGM due to a UV background. The number of de- 
tectable fossil H II regions would then equal the number 
in equation (30) multiplied by the inverse of the quasar 
duty cycle ~ 1% (Wyithe & Loeb 2003b), or about 200 
per field. 

LOFAR should also be sensitive to H II regions around 

less luminoiis quasars. The brightness temperature of the 
emission ring, or of the IGM outside the H II region has a 
mild dependence on the quasar luminosity, and its diame- 
ter scales only with the one third power of the luminosity. 
H II regions around less luminous quasars should be more 
numerous by a factor proportional to L'^"'"^. For example, 
detection of redshifted 21cm emission around quasars with 
luminosities a factor of 10 smaller than those of the bright 
SDSS quasars require a resolution that is better by a fac- 
tor of 10^/^, corresponding to an integration time that is 
longer by a factor of lO^/"* ^ 4.6. However, the number of 
active high redshift quasars per field powering such emis- 
sion rings is likely to be of order 10^ ~ 200, while the 
number of fossils regions could be ^ 2 x 10"*. 

Finally, we briefly discuss the example of a redshifted 
21cm signature around a high redshift quasar presented in 
5b of Tozzi et al. (2000). Their example shows a ring of 
redshifted 21cm emission inside a ring of redshifted 21cm 
absorption. This morphology is unlikely to be observed 
because the observation of an absorption ring requires the 
absence of a UV background, and the provision of coupling 
of the spin and kinetic temperatures through Lya photons 
from the quasar itself. As we have shown, Lya photons 
from the quasar are unable to couple the spin and kinetic 
temperatures. An emission ring enclosed within an IGM 
that absorbs everywhere, as shown in figure 5 is the cor- 
rect signature of a Stromgren sphere expanding into a cold 
IGM (provided that there is a UV background). 

6. STROMGREN SPHERES AND THE GLOBAL SIGNATURE 
OF REIONIZATION 

A very exciting possibility for instruments like LOFAR 
is to search for the emission (or absorption) of redshifted 
21cm radiation from the neutral IGM (averaged over large 



scales) before reionization. By finding the frequency at 
which this emission (absorption) cuts off, one might deter- 
mine the redshift of reionization (e.g. Shaver et al. 1999). 
If reionization took place over an extended period of time 
as suggested by recent work (Barkana & Loeb 2003b), then 
this global reionization signature will be difficult to detect 
for two primary reasons: (i) because it will be superim- 
posed on a smooth but varying foreground spectrum, and 
(ii) calibration needs to be performed over a frequency 
range that is broader than tlic^ detector band-pass. 

Arc-minute resolution observations of quasar Stromgren 
spheres would provide an alternative method to search for 
a step in the brightness temperature of the IGM due to 
reionization. This is because the inner part of the H II 
region provides a calibration of the brightness tempera- 
ture of the radio foreground that is unavailable along other 
lines of sight. Note that 21cm emission from the quasar 
itself does not affect this calibration, since it effectively 
just adds to the foreground. Observations along the line- 
of-sight to the quasar would include emission within the 
spectral ranges corresponding to both the ionized region 
and the neutral gas within the same band-pass. Redshifted 
21cm observations of quasars at a series of redshifts would 
therefore allow determination of the brightness tempera- 
ture of the IGM at different cosmic times, and a detection 
of the onset of reionization. Redshifted 21cm observation 
of highest-redshift SDSS quasars will immediately deter- 
mine whether the IGM was significantly neutral at ^; > 6.3. 

7. DISCUSSION 

Observations of the reionization epoch in redshifted 
21cm radiation with proposed instruments like LOFAR 
will offer a unique opportunity to study the state of the 
IGM at the end of the dark ages. The warm neutral hy- 
drogen outside the H II regions surrounding bright high 
redshift quasars is expected to emit at 21cm. The bright- 
ness temperature contrast of this emission will probe the 
neutral fraction (with a lower contrast implying a lower 
neutral fraction), although the inference will be subject 
to uncertainties in the temperature of the warm IGM. Al- 
ternatively, the relativistic expansion of the H II regions 
provides a geometric method to determine the neutral frac- 
tion as well as the quasar lifetime. 

We have combined the observed Lya absorption signal 
with expectations for redshifted 21cm signatures of the 
H II regions around the highest redshift SDSS quasars. Be- 
cause the H II regions around the known z > 6.3 quasars 
are at a stage of relativistic expansion, the line-of-sight 
and transverse radii (measured through Lya absorption 
and redshifted 21cm emission, respectively), correspond to 
two different epochs in the evolution of the H II regions. 
We have found that measurements of these two radii can 
be used to determine the neutral fraction of the IGM with 
only a mild dependence on the anisotropy of the quasar 
emission (see Fig. 2). Assuming isotropic quasar emission 
and the absence of a fossil H II region prior to the quasar 
activity, the observation of a redshifted 21cm emission ring 
with a radius greater than 2.5Mpc (0.5Mpc) around the 
known z > 6.3 quasars, would provide a lower limit on 
the neutral fraction of a;Hi > 0.3 (xm > 0.05). Emission 
concentrated close to the line-of-sight and the existence of 
fossil H II regions around the quasar host galaxy would 
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imply even larger neutral fractions. 

At low redshifts, quasar lifetimes are bracketed by a 
variety of techniques to the range ^ 10^-10*yr (e.g. Yu 
& Tremaine 2002; Martini & Weinberg 2001; Haiman & 
Hui 2001; Jakobscn, Janscn, Wagner & Rcimcrs 2003), 
with preferred values near 10'' years (e.g. Martini 2003). 
The lifetime required for the H II region to expand to 

— 1/3 

its observed size is proportional to x^^ . The geomet- 
ric determination of the neutral fraction can therefore be 
directly used to determine also the quasar lifetime mea- 
sured from the combination of Rp and Xm- Thus, Lya 
and redshifted 21cm measurements of the expanding H II 
regions around the highest redshift quasars will allow a di- 
rect determination of quasar ages on a case by case basis 
and constrain quasar formation models (e.g. Haiman & 
Loeb 2001). 

Based on current number counts of bright high redshift 
quasars and the preliminary design of LOFAR, we fore- 
cast that each LOFAR field (of radius ~ 11°) will contain 
~ 2 active quasars as bright as the observed SDSS quasars 
within one band-pass of frequency (^ 4MHz) . An integra- 
tion time of ^ 100 hours should be sufficient to detect the 
rings of redshifted 21cm emission from the warm IGM out- 
side the quasar H II regions if xhi = 1. In the presence of a 
strong cosmic Lya background, these rings should remain 
visible after the quasar has turned off. This would result 
in the detection of a large number of fossil H II regions 
that exceeds the number of active quasars by the inverse 
of the quasar duty cycle (~ 100). 



Arc-minute observations of quasar Stromgren spheres 
will provide an alternative method for measuring the spec- 
tral step that makes the so-called global signature of reion- 
ization (Shaver et al. 1999). The brightness temperature 
measured in the spectral range corresponding to the ion- 
ized gas within the H II region provides a calibration of 
the foreground emission at the same redshift. This cali- 
bration will help overcome the difficulties associated with 
measuring the global signature of an extended reioniza- 
tion epoch (as suggested by the large electron scattering 
optical depth measured by WMAP; Kogut et al. 2003) on 
top of a frequency dependent foreground. Detection of the 
signal for a statistical sample of quasars would help cali- 
brate the expected scatter in the reionization redshift of 
different regions in the universe (Barkana & Loeb 2003b). 

Finally, we emphasize that the redshifted 21cm observa- 
tion of existing SDSS quasars will determine whether the 
IGM was significantly neutral at z > 6.3 as suggested by 
their observed Lya spectra (Wyithe & Loeb 2004). When 
combined with the high value of the optical depth for elec- 
tron scattering measured by WMAP (Kogut et al. 2003), a 
significant neutral fraction at 2; ~ 6.3 would imply a com- 
plex reionization history and possibly an early reionization 
phase driven by Population-Ill stars. 
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